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SUMMAEY 

An approximate method, "based on a simplified form of the con- 
tinuity relation, is developed to predict the location of detached 
shock waves ahead of two-dimensional and axially symmetric bodies. 
In order to reduce the problem to an equivalent one -dimensional 
form, it was assumed that: (1) The form of the shock between its 

foremost point and its sonic point is adequately represented by 
an hyperbola asymptotic to the free -stream Mach lines; and (2) the 
sonic line between the shock and the body is straight and inclined 
at an angle that depends only on the free -stream Mach number. With 
these assumptions, the location of the shock relative to the body 
sonic point Is independent of the form of the nose or leading edge 
ahead of the sonic point and becomes a single-valued function of 
the Mach number. A simple geometric method for estimating shock 
location is also presented, but this method agrees less closely 
with experiment than the continuity method. 

When the location of the detached wave and the sonic line 
relative to the body are known, the drag of the portion of the 
body upstream of its sonic point can be estimated from the momen- 
tum change of the fluid that crosses the sonic line. Comparison 
with available experimental results Indicates that the estimated 
drag coefficients are good approximations except for very blunt 
bodies and for free -stream Mach numbers close to 1.0. 

The simplified continuity method Is also applied to unsym- 
netric two-dimensional, bodies at angle of attack and to two- 
dimensional and axially symmetric supersonic inlets with detached 
waves due to spillage over the cowl. For supersonic inlets, the 
relation between spillage and shock location is presented and the 
additive drag due to spillage is estimated. 
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INTRODUCTION 

When a region of subsonic flow occurs in a supersonic flow 
field, many simplifying assumptions are generally required to make 
theoretical analysis feasible. Without such assumptions, the flow 
can be constructed only by means of lengthy numerical methods 
based on differential equations of fluid mechanics. Such compu- 
tations differ for each configuration investigated; general trends 
and important parameters are therefore not easily discernible. 

A reasonable set of assumptions for simplifying the analysis 
of supersonic flow with detached shock waves should be consistent 
with the following characteristics of such flow fields: The entropy 

varies from the free -stream value to the value behind a normal 
shock, and the velocity varies from the free-stream value to zero 
at the stagnation point. Thus, if the free-stream Mach number is 
considerably greater than unity and if the subsonic region is a 
significant part of the entire disturbed flow field, analyses based 
on such assumptions as irrotational flow, small differences between 
local and free-stream velocity, or incompressible flow in the subsonic 
region can no longer be expected to yield valid results. For such 
problemd, other simplifying assumptions must therefore be formulated. 

An analysis based on assumptions concerning the form of the 
boundaries of the subsonic region rather than on the nature of the 
flow variables is presented herein. In this analysis, which was 
made at the NACA Lewis laboratory, both plane and axially symmet- 
tric bodies were considered. The form of the detached wave is 
assumed to be only secondarily influenced by the form of the body 
ahead of the sonic point, and the sonic line between the shock and 
the body is assumed to be straight and inclined at an angle that 
depends only on the free-stream Mach number. The continuity rela- 
tion is then applied to this simplified picture to obtain the loca- 
tion of the shock relative to the sonic point for plane and axially 
symmetric bodies. The drag is estimated from the momentum change 
of the fluid that passes the sonic line. The method is also 
applied to plane unsymmetric bodies at angle of attack and to 
two-dimensional and axially symmetric supersonic inlets with 
spillage over the cowl. Results are compared, where possible, 
with available experimental data. 
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SYMBOLS 


The following symbols are used, some of which are illustrated 
in figure 1: 

A area 

a speed of sound 


B **’ 

b 

C 


k = 


L 


M 

P 

P 

V 

x 


y 

P 

7 



x- coordinate of foremost point of body 


P^P tan <p s - aJ p 2 tan 2 cp g - lj 


drag coefficient of portion of body ahead of theoretical sonic 
point (based on area indicated by subscript) 


/ 8 \ O 

1.2679 — ) - 1 

VoA Po 


tm 


0 


distance between vertex of detached shock wave and x, 


SB 


free -stream Mach number 
stagnation pressure 
static pressure 
local velocity 

coordinate in stream direction 

distance from foremost point of detached shock to intercept 
of its asymptote on x-axis 

coordinate perpendicular to stream direction 

-(/m 2 - 1 

ratio of specific heats (1.4 for air) 
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8 local angle of "body contour relative to x-axis 

Tj angle between sonic line and normal to free -stream direction 

9 half -angle of wedge or cone 

cone half -angle for which shock becomes detached 
A angle of streamline relative to x-axis 

wedge half-angle for which shock becomes detached 
p density 

CT isentropic contraction ratio from free stream to sonic velocity 

T fraction of maximum possible inlet mass flow that passes 
outside cowl 

cp local inclination of detached shock relative to x-axis 
Subscripts : 

0 free -stream conditions 

1,2 upper and lower portion of unsymmetric body contour, 
respectively 

c centroid of stream tube passing sonic line 

cr critical conditions 

S sonic point of detached shock 

SB sonic point of body 

s conditions along sonic line 


ANALYSIS 

The representation of flow with detached shock waves to be used 
in the present analysis is shown in figure 1. The region of subsonic 
flow Is bounded by the portion of the detached wave below its sonic 
point S, by the portion of the body contour upstream of its sonic 
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point SB, and by the sonic curve between S and SB. If the form 
of the detached wave is known, the location of S on the shock can 
immediately be determined, because for a given Mq the angle qpg 

for which sonic velocity exists behind the shock is known from shock 
theory. The flow direction Xg at this point Is then also known. 

The sonic curve is shown in figure 1 as a straight line. This 
approximation is used throughout the analysis, although more exact 
computations, as well as experimental results, indicate that the 
form may depend considerably on the shape of the nose or leading 
edge (references 1 and 2). To the degree of approximation of the 
present method, however, such variations appear to be unimportant. 
With the simplified picture shown in figure 1, approximate expres- 
sions can be derived for the shock location relative to the body 
sonic point and for the drag of the portion of the body in the 
subsonic-flow region. 


Location of Body Sonic Point 

Evidence is available to show that, for bodies with sharp or 
well-defined shoulders (fig. 2(a)), the sonic point is located at 
the shoulder (reference 2). For more gradually curved bodies, such 
as ogives, the location of the sonic point can be estimated by 
an extension of a method suggested by Busemann (reference 3). This 
method locates the shoulder at the point where the contour of the 
body is inclined at the wedge angle or cone angle corresponding to 
shock detachment (fig. 2(b)). If the sonic line is approximated by 
a straight line, the location of the shoulder and the sonic point 
coincide. These results are analogous, in some ways, with the 
flow past the throat section of a supersonic nozzle, for which the 
location of the mw.ri mm constriction (shoulder) and the sonic point 
coincide if the throat is sharp or if the flow is treated as one 
dimensional. 

The choice of the point of tangency of the body with a line 
inclined at the detachment angle aB the approximate location of the 
sonic point has some theoretical Justification for two-dimensional 
flow. W. Perl has pointed out to the author that for flow near 
sonic velocity and for small perturbation velocities, the deriva- 
tive dX/dV taken along the normal to a streamline, vanishes to a 
first approximation at the intersection of the sonic line with the 
shock wave. This result implies that X is relatively constant 
along the sonic line and hence, that the inclination of the body 
at its sonic point is approximately equal to Xg. Inasmuch as X g 
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la slightly smaller than A d , the sonic point would he located 

slightly downstream of the shoulder of the body; whereas the dis- 
cussion given in reference 3 shows that the sonic point lies ahead 
of the shoulder. These differences. However, are insignificant for 
the present analysis, because the sonic point and the shoulder are 
assumed to coincide. 

Although the theoretical basis for locating the sonic point 
by means of the detachment angle is far from complete, particularly 
for axially symmetric bodies, the procedure appears reasonable if, 
as assumed in the present analysis, the form of the nose or lead ing 
edge ahead of the shoulder (or sonic point) has only a secondary 
effect on the form and the location of the sonic line. This sonic line 
should therefore have almost the same form for all bodies as, for 
example, for a wedge or a cone with half -angle only slightly greater 
than the detachment angle. Because the sonic point for these bodies 
is known to be located at the point of tangency with the detachment 
angle, the procedure should be valid for all bodies at the same 
Mach number within the degree of approximation of the present method. 
This procedure will therefore be used to determine body sonic points 
for comparison of theoretical and experimental results. 

The reasoning used in the preceding paragraph leads to a very 
simple geometric method for predicting the location of detached 
shock waves. A wedge or a cone with included angle large enough to 
cause shock detachment and with well-defined shoulders is considered 
(fig. 2(a)). If b-x Q is the distance from the vertex of the cone 

or wedge to the foremost point of the detached wave, then 


b-x c 

y SB 


L 

y SB 


cot 6 


( 1 ) 


If 6 is now set equal to the detachment angle, then b-XQ becomes 
zero, and the expressions for shock location become: 

for plane flow. 


-h— = cot Aj 
y SB 

and for axially symmetric flow, 


( 2 ) 


= cot e d 


L 

y SB 


(3) 
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Because, by the present assumptions, L/ygB is constant for a given 
Mq and independent of the form of the nose, equations (2) and (3) 

represent estimates of the location of the detached shock wave for 
all bodies at zero angle of attack. This estimate should be par- 
ticularly valid for bodies only slightly blunter than the wedge or 
cone for which detachment occurs. The following derivation, how- 
ever, which is based on the continuity relation, should provide a 
better estimate of the average shock location for more general 
blunt bodies. 


Location of Detached Shock Wave 
(Continuity Method) 

Assumed form of detached shock wave. - The most typical char- 
acteristics of detached waves are that: (l) They are normal to the 
free stream at their foremost point; and (2) they are asymptotic to 
the free-stream Mach lines at large distances from their foremost 
point. A simple curve that has these characteristics is an hyperbola 
represented by 


Py = ^x 2 -x Q Z (4) 

where 0 is the cotangent of the Mach angle and x Q is the dis- 
tance from the vertex of the wave to the intersection of its 
asymptotes. For the purpose of the analysis, the hypothesis is made 
that to the degree of approximation required, all detached waves in 
the region between the axis and the sonic point S may be 
represented by equation (4). 

With this form of detached wave, the angle between the stream 
direction and the tangent to the shock at any point is obtained 
from (fig. l) 


= tan cp 
dx 


4 


.A) 


x 2 -x 2 

X Xq 


x Q 2 +0 2 y 2 

0 2 y 


(5) 
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The location of S is then 



If the y-coordinate of the body sonic point is used as the refer- 
ence dimension, then the form of the shock wave Is given by 
(equation (4)) 



where, from equation (6), 



From equation (6a) the dimensionless location of the shock sonic 
point is 



The distance from the foremost point of the shock to the x-coordinate 
of the body sonic point is 

L _ ^SB X 0 

y SB y SB y SB 


( 10 ) 
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where, from figure 1, 

fsB _ fs_ + / ys 
ysB • ysB \ys® 

If equations (8), (9), and (ll) are combined to eliminate all 
unknown coordinates except yg and L, equation (10) becomes 

y 

JL = JL (C + tan n) - tan n (12) 

y SB y SB 

where 


C = P (^3 tan cp s - /\jp 2 tan 2 cp g -lj (13) 

Inasmuch as P and cpg are known for any given free-stream 

Mach number, only the quantities y^/^SB *1 remain to be 

determined to predict the relation between the sonic point on a 
body and the location of its detached wave. 

Applicati on of continu ity equ a tion. - In order to determine 
the quantity yg/yggj the continuity relation is applied to the 

fluid that passes the sonic line. The integral form of this rela- 
tion is useless for the present method, inasmuch as the distribu- 
tion of the flow variables along the sonic line is unknown. The 
distribution of stagnation pressure immediately behind the shock, 
however, is known. An appropriate average value of this quantity 
is that existing along the streamline which represents the mass 
centroid of the fluid passing the sonic line. This centroid 
streamline enters the shock wave at y c =yg/2 for plane flow and at 

y c =2yg/3 for axially •syrnmetrj c flow. The shock angle correspond- 
ing to these values of y can be obtained from equations (5) and 
(6). Because the stagnation pressure remains constant along each 
streamline behind the shock, the value of P 8>c will remain 

unchanged between the shock wave and the sonic line. Because the 
total temperature is constant, the simplified continuity equation 
may be written as 


1) tan n (11) 

/ 
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A 0_ Wo _/ P 9 \ < pV >cr 
A s P 0 V 0 \ p o/ c _ P 0 V 0 _ 



(14) 


where a is the contraction ratio required to decelerate the free 
stream to sonic velocity isentropically. In terms of the coordi- 
nates of the sonic points, equation (14) becomes, for plane flow 


or 


y S~ y SB 
cos T) 



b^s 


— = (l - B cos T}) 

ysB 


(15) 


and for axially symmetric flow 


y S 2 ~ y SB 2 
cos n 




or 


£s 

y SB 


(1 - B cos q) 


1 

2 


(16) 


The appropriate values of tj to be used in each case remain 
to be established. On the basis of the previously cited analogy 
between one-dimensional channel flow and flow with detached shock 
waves, the sonic line is assumed to be normal to the average flow 
direction in its vicinity. At S the inclination of the flow is 
known to be Ag both for plane and axially symmetric flow; 

whereas at SB the inclination is assumed to be A a for plane 

flow and 0^ for axially symmetric flow. If the arithmetic 

mean of the inclinations at the two extremities is used, the 
appropriate expressions for t] become: 
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for plane flow, 

1 - 5 (\l ♦ Vi) (17) 

and for axially symmetric flow, 

1 = \ (6 d + V (18) 

Because Ag differs only slightly from A^, the inclination of 

the sonic line for plane flow will be assumed to be simply tj = 

Values of Js^SB stained from equations (15) and (16) are shown 

as functions of Mq in figure 3. The variation of Xg/Vgg with 

Mq, as computed from equation (8), is also shown. The values of 

Ag and required to determine yg/y^ were obtained from the 

shock charts of reference 4 and are plotted against Mq in figure 4. 

The shock angle at the sonic point cpg and the shock angles at the 

mass centroid for two-dimensional and axially symmetric flow are also 
shown in figure 4. In figure 5, the resulting total-pressure ratios 
(P s /P q) 0 are given together with the values of a and C used in 

equations (13) through (16). 


Comparison of Theoretical and Experimental Shock Location 

A comparison of the shock form and location estimated by the 
continuity method with the experimental results of reference 2 is 
shown in figure 6. The data of reference 2 were obtained with a 
free jet with an outlet Mach number of 1.7. The flow field in the 
vicinity of the models was reconstructed from interf erograms . For 
comparison of theory and experiment, the foremost point of the 
detached wave was used as the common point. The theoretical con- 
figuration is independent of the body form and is therefore the 
same for each of the three bodies shown. For the sphere and the 
cone, experimental sonic points are given in reference 2 and the 
y-coordinate of these points was used as the reference dimension. 
For the projectile, however, the actual sonic point was not deter- 
mined in reference 2 and the theoretically estimated location of 
this point was used as the reference dimension. The assumed shock 
form and the estimated distance from the vertex of the shock to the 
body sonic point are shown to be approximately in agreement with 



12 


NACA TN 1921 


experimental results for all three bodies. Increasing bluntness of 
the nose appears to shift the upper portion of the detached wave 
away from the body, but in view of the wide range of nose forms 
this effect appears to be small. The shock form for the projectile 
may be somewhat in error in figure 6 due to the relatively small 
scale of the subsonic portion of the diagram presented for this 
body (reference 2). 

The theoretical variation of L/ygg with Mach number is shown 

in figure 7 for plane and axially symmetric bodies, together with 
available experimental values from references 1, 5, and 6, and the 
numerically computed values of reference 1. The data obtained from 
reference 6 were converted to the present parameters from a small 
sketch of the model and may therefore be somewhat in error. The 
values obtained by the continuity method are seen to agree more 
closely, in general, with experimental results than the values 
obtained by the geometric method. As previously stated, however, 
the geometric method should yield more accurate results for bodies 
only slightly blunter than the wedge or cone for which shock detach- 
ment occurs. This prediction is bom out to some extent for the 
two-dimensional bodies, where the value obtained in reference 1 for 
the wedge agrees closely with the value obtained by the geometric 
method, and the continuity method appears to average the results 
obtained for the wedge and the flat -nosed body. For the axially 
symmetric bodies no such comparison is available, although fig- 
ure 6 shows that L/ygg tends to increase as the body becomes less 

blunt. These results indicate that the continuity method leads to 
an average value of L/ygg, whereas the geometric method yields 

the maximum value of this quantity for a given Mq. 


Drag Upstream of SB 

Because the form and the location of the shock wave and the 
sonic line are assumed to be independent of the shape of the nose 
or the leading edge, it follows that the drag coefficient of the 
portion of the body upstream of its Bonic point is, by the present 
method, also independent of the body form. This portion of the 
drag can be determined from the momentum theorem. If the stream- 
lines are assumed to be normal to the sonic line and if average 
values of the flow variables at the sonic line are again used, the 
momentum equation may be written as (fig. l) 
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OSB 


(p-p 0 )dA sin 6 - P 0 V 0 2 A 0 -(PgV 8 2 ) c A 8 cos 1-(P B#0 -P 0 )A a cos 1 


from which the drag coefficient based on Aq becomes 
r SB 


< C d) 


(p-pJdA sin 5 
lb 0 

— P V 2 A 
2 0 0 A 0 


(19) 


r 

<i - 

( v„\ + 

fts) - i 

V p oA 

B COS T1 

_ 2 + 

2 


V 

P V 

L oo 

* 

o 

1 

0 

y 


2 1 + 

0 


1 - (7+1) 


(l sV 

w. 


(20) 


Using the value of A 8 /Aq frost equation (14), equation (20) 

becomes 


r 


(cj » 2 ( l + 2-222JL 


> P ° M 2 I 

7 ^ M o L 


p o P 0 


Ps 


p 0 yl ■ p s 
' c 


(7+1) 


( 21 ) 


where ^ = 0.5283. 

The drag coefficient based on the area of the body at its 
sonic point is : 
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for plane flow, 


(0 ->\. ' ^ &) 
and for axially symmetric flow, 



( 22 ) 


(23) 


If the appropriate values of tj, (P S /P Q ) , and y s /y^ are used, 

the expressions for the drag coefficient "become idential for plane 
and axially symmetric flow: 


where 



' 1 - kB cos q 
il - B cos tj j 



(24) 


(25) 


These drag coefficients are plotted against Mq in figure 8, 

together with the experimental values obtained from reference 2 by 
numerical integration of the pressure coefficients given therein 
and the theoretical values for flat-nosed bodies and for a wedge 
computed in reference 1. The theoretical values obtained by the 
present method are not expected to be valid near Mq = 1.0, because 

the subsonic field becomes very large relative to the body 
( y s/ y SB — * 09 ^ en< l small variations in the assumed form and 

inclination of the sonic line have a large effect on the computed 
drag. Equation (24) becomes indeterminate for Mq = 1 because 

cos tj, B, and k are all unity for sonic free-stream velocity. 

The agreement with experimental data for the cone and the sphere 
at Mq = 1.7 is seen to be good. The drag coefficients computed 

for the flat-nosed bodies in reference 1, however, fall considerably 
above the values obtained by the present method. Part of this 
disagreement is probably due to the somewhat extreme form of the 
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sonic line assumed for these bodies In the computations of refer- 
ence 1. This assumed sonic line meets the face of the body per- 
pendicularly at the shoulder, whereas the actual sonic line is 
probably concave downstream. (See fig. 6.) Some disagreement 
should probably be expected, however, Inasmuch as the momentum 
integral used to obtain the drag coefficient by the present method 
is more sensitive to differences in the stream angle near the 
sonic line than the continuity equation used to predict shock 
location. The present procedure is therefore most valid when the 
streamlines are expected to be almost normal to the sonic line and 
when the sonic line is not expected to have large curvature. 

The maximum possible drag coefficient, obtained by assuming 
that the stagnation pressure behind a normal shock acts on the 
entire nose, is plotted in figure 8 for comparison. 


Unsymmetric Bodies and Effect of Angle of Attack 
for Two-Dimensional Flow 

If the method used for determining the location of SB for 
symmetric bodies is applied to unsymmetric bodies or to bodies 
at angle of attack, the magnitude of the coordinates of the upper 
and lower sonic points will, in general, differ (fig. 9). It is 
apparent that portions of figure 9 above and below the center line 
between the two tangent lines are each similar to the configuration 
shown in figure 2(b). Such a representation, however, can be con- 
sidered valid only for bodies that can be completely described in 
two dimensions. If the reasoning used to locate the sonic point 
for symmetric bodies at zero incidence Is followed, then varia- 
tions of the portion of the body between SB,1 and SB, 2 will not 
appreciably alter the location of the sonic points or the form and 
the location of the detached wave and the sonic line. The conti- 
nuity method can be applied separately to the portions of the flow 
field above and below the zero streamline (that is, the streamline 
that reaches the stagnation point on the body) . The location of 
this streamline need be known only in the free stream ahead of the 
detached wave. The assumption that this streamline must pass 
through the intersection of the two tangent lines (fig. 9) is con- 
sistent with the reasoning used to locate the sonic points. Because 
the upper and lower portions can be treated separately, equations (15) 
and (16) lead to the result 


y S,l = y S,2 = £s_ 
y SB,l y SB,2 y SB 


( 26 ) 
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Therefore, from equations (8) and (12), 



y SB,l y SB,2 y SB 


and 

L . L B L 
y SB, 1 y SB, 2 y SB 


(27) 


(28) 


The equations for the portions of the shock wave above and 
below the x-aiis are: 


Upper 


Lover 



(29) 


(29a) 


An example of the resulting configuration for an unsymmetric 
two-dimensional body is shown in figure 10 for y SB,2 = 0 ' 5y SB . 1 . 

The values of and * 0 / y SB were obtained from figure 3 and 

L/ygg was obtained from the results of the continuity method 

(fig. 7). A slight discrepancy appears in the location of the 
vertex of the detached wave due to the separate construction of the 
upper and lower portions of the configuration, but this discrepancy 
is well within the accuracy of the approximations. Use of the 
geometric method would locate the origin of the shock at the point 
of intersection of the tangent lines from the sonic points. This 
method would avoid the discontinuity at the vertex of the shock, 
but would fail to satisfy the continuity equation. 

From equations (21) and (22), the drag coefficient of the lower 
and upper portions of the body are seen to be dependent, for a 
given Mq, only on y g /y gg, so that the drag coefficient for the 
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portion of the unsymmetric "body ahead of the sonic points based 
on ( y SB 1 + ^SB 2 ) identical to the drag coefficient based on 
Agg for the symmetric body at zero angle of attack. This result 

is also evident from the consideration that the momentum change of 
the fluid that passes the sonic line depends only on the Mach number 
of the free stream. The total mass flow past the sonic lines 
depends only on Agg, but the portions of this mass flow above and 

below the zero streamline are not equal for unsymmetric bodies or 
for bodies at angle of attack. 

In addition to the limitations mentioned in connection with 
the estimated drag coefficients at zero angle of attack, the fact 
that regions of negative pressure coefficient may occur when the 
stagnation streamline fails to intersect the foremost point of the 
body must also be considered. The drag coefficients estimated by 
the present method correspond to the integral of the pressure forces 
over the zero streamline and are independent of the shape of this 
streamline ahead of the sonic point. Local separation regions near 
the foremost edge, however, may alter the subsonic portion of the 
flow sufficiently to i n validate estimates of drag based on a pre- 
scribed form of the sonic line and an assumed flow direction. The 
procedure for estimating the location of the shock, however, should 
not be greatly in error, because the continuity relation used to 
predict this location is relatively insensitive to variations in 
the form of the sonic line or the inclination of the velocity vector. 


Shock Location and Additive Drag for Supersonic 
Inlets with Spillage over Cowl 

The method of the preceding sections may be easily extended 
to the problem of estimating the shock location and the additive 
drag as functions of the fraction of the maximum mass flow that is 
spilled over the cowl of two-dimensional or axially symmetric 
inlets. The simplified picture corresponding to this problem is 
shown in figure 11, where the y-coordinate is now measured from the 
axis pf symmetry of the inlet. The coordinate y m denotes the 

free-stream location of the streamline that separates the mass 
entering the inlet from that passing outside the inlet. The inter- 
section of this streamline with the detached shock wave is assumed 
to be the origin of the hyperbolic portion of that shock wave. In 
terms of the coordinates shown in figure 11, the spillage t Is 
defined by the following expressions: 
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Two-dimensional inlets 


T ■ 1 - (30) 

y SB 

Axially symmetric inlets 



The drag due to spillage is defined as the integral of (p-Pg) 

along the streamline that bounds the entering fluid. This drag is 
similar to the additive drag defined by Ferri for nose inlets with 
projecting central bodies, for which the streamlines may be deflected 
a considerable distance ahead of the cowl. 

The sonic point SB can be estimated by the procedure used 
for closed bodies. For most cases of interest, the lip of the cowl 
will be sharp or very thin relative to the inlet cross section, so 
that SB will be very close to the foremost point of the lip. Here 
the flow rapidly expands to supersonic velocities. Separation 
regions resulting from this expansion are again neglected and may 
considerably reduce the additive drag due to spillage. The relation 
to be derived between spillage and shock location, however, should 
not be greatly affected by local separation regions. 

Other complications usually encountered with supersonic inlets 
have also been ignored. Oblique shocks from a protruding central 
body, for example, are assumed to have negligible effect on the 
relation between mass-flow spillage and drag and on the form of the 
shock wave above y m . These simplifications should yield valid 

estimates when open-nosed bodies without central bodies are con- 
sidered or when the effect of the central body is felt only below 
the origin of the hyperbolio portion of the detached shock. The 
analysis is independent of the nature of the shock below this 
origin. 

The hyperbolic portion of the shock wave is now represented by 
the equation 


JSB 




( 31 ) 
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where 



0 2 tan 2 cPa-1 



(32) 


Equations (9), (10), and (ll) for the x-coordlnates of the sonic 
points and for the shock location remain unchanged by translation 
of the origin. When equation (32) is used in place of equation (8) 
the expression for the shock location (equation (10)) becomes 


-JL. = (C + tan tj) - C - tan (33) 

y SB y SB y SB 


Two-dimensional inlets. - For two-dimensional flow, the ratio 
of the sonic-line flow area to the free-stream area of the stream 
tube is (fig. 11) 


^s _ y S-ySB = Ao \ = B 
A 0 ’ ( y s- y m> 008 A s V p 84 


(34) 


so that 



1 - 


B cos A 

y SB 


1 - B cos X 

s 


(35) 


If this value is substituted in equation (33), the expression for 
the location of the vertex of the detached wave becomes 


L 

y SB 



Jm\ /C + B sin X g \ 

W v 1 ' B cos X S / 


(36) 


When no air passes through the inlet (y m =0), this expression is 
found to reduce to the value obtained for closed symmetric bodies . 


The drag coefficient based on Aq (equation (21)) is inde- 
pendent of the coordinates, so that the drag coefficient for the 
inlet is obtained by multiplying (C D ) by the ratio Aq/A^, 

A 0 


and becomes 
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(Cd) Asb " (Cd) A 0 ( y SB ’ 7 Sb) “ (Cd) A 0 Cl " 1 + T ) (37) 

Equations (35) and (36) indicate that the distance of the shock 
from the inlet lip and the drag coefficient for a given Mq 

increase linearly with the amount of air that passes outside the 
inlet . 


Axially symmetric Inlets. - The ratio of sonic area to free- 
Btream area for axially symmetric inlets becomes 


ys 2 -ysB 2 


0 (y s 2 " y m ' cos ^ 


= B 


so that 


(38) 


£s_ 

7 sb 

For axially symmetric inlets, however, the value of ^ to be 
used may be expected to vary with the value of y m /ycjg. When this 

ratio is close to unity, phenomena near the cowl lip should be 
almost two dimensional j whereas for large spillage the flow becomes 
similar to that obtained with closed bodies of revolution. Inas- 
much as the cases of relatively small spillage are probably of 
greater interest, it will be assumed that the flow near the cowl 
lip is two dimensional so that q = Xg. Similarly, the two- 

dimensional value of total-pressure ratio at the mass centroid will 
be used. By substituting equation (39) into equation (33), the 
shock location is found to be 


1 ' fl J B 003 11 


1 - B COS TJ 


(39) 



Ji) B 008 *s 


1-B cos X c 


(C + tan Ag) 



— tan 


( 40 ) 
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ratio 
area, 

+ T 

(41) 

where q in equation (21) is replaced by for y m ~ ygg. 

The variations of shock location and drag coefficient with T 
are shown in figure 12 for two-dimensional and axially symmetric 
inlets at Mq = 2.0. For axially symmetric inlets, the values 

obtained using q = - (0, + X g ) are also plotted for comparison. 

2 

The value assumed for q evidently has little effect on the esti- 
mated drag or the shock location. 


Equation (21) for (C D ) may again be multiplied by the 

A 0 

Aq/A<-jb "to obtain the drag coefficient based on the inlet 
Thus, 


(0 d ) Asb = (C D ) J 


— 

V'SB, 


( £ m 

VSB> 


» (c D ), 


*0 


*L) - 1 

V y SB/ 


SUMMARY OF METHOD 

An approximate method has been developed for predicting the 
location of detached shock waves ahead of plane or axially 
symmetric bodies and for estimating the drag of the portion of 
the body upstream of the theoretical sonic point. The method is 
based on the continuity relation, which is applied to the air that 
passes the sonic line. The main assumptions were that: (l) The 

form of the detached wave from the axis to its sonic point is ade- 
quately represented by an hyperbola asymptotic to the free-stream 
Mach lines; and (2) that the sonic curve between the shock and the 
body is a straight line. These assumptions imply that the form and 
the location of the shock and the sonic line are only secondarily 
influenced by the form of the body ahead of its sonic point and that 
the shock location relative to the body sonic point is primarily a 
function only of the free-stream Mach number. The drag of the body 
to the sonic point was estimated from the change of momentum of the 
air that passes the sonic line and was consequently also independent 
of the form of the nose or the leading edge in the present approxi- 
mation. 
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Comparison with experimental results indicates that the assump- 
tions made are valid to good approximation for predicting the form 
and the location of detached waves, but may be oversimplifications 
when it is desired to predict the drag of extremely blunt bodies. 

The method also fails in the vicinity of sonic free-stream velocity. 
The effects of angle of attack and unsymmetric body contour are 
discussed for two-dimensional flow, but no experimental results are 
available to check the theoretical predictions. 

Extension of the method to the problem of estimating the addi- 
tive drag and the shock location for supersonic inlets with spillage 
over the cowl indicated that the drag, both for two-dimensional and 
axially symmetric inlets, increases linearly with the percentage 
of the m*Tl mmn inlet mass that passes outside the cowl. The dis- 
tance between the shock and the inlet lip also Increases linearly 
with this percentage for two-dimensional inlets and almost linearly 
for axially symmetric inlets. 


Lewis Flight Propulsion Laboratory, 

National Advisory. Committee for Aeronautics, 
Cleveland, Ohio, June 1, 1949. 
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Figure 1. - Representation of flow with detached shock wave 
and notation used In analysis. 
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(a) Well-defined shoulders (plane or axially symmetric). 






Plane bodies 


Axially symmetric bodies 


(b) Gradually-curved bodies. 


Figure 2. - Approximate location of body sonic point. 
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Figure 4. - Variation with Mach number of several angles used in analysis^ 
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Figure 5. - Variation with Mach number of acne parameters used In analysis 
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Figure 7. - Theoretical and experimental shock location. 
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Drag of portion of body upstream of theoretical sonic point 
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Figure 9. - 


Determination of sonic points for plane unsymmetric 
bodies at angle of attack. 
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Figure 10. - Predicted relation between detached shock body 
sonic points for two-dimensional body when y m o.5v 1 
Mq * 2.0. SB,2 ° oB,l 
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Figure 11. - Simplified picture used for analysis of flow past 
supersonic inlets with spillage. 



Figure 12. - Shock location and drag as function of spillage* Mq s 2JX 
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